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AN EXPERIMENTAL INVESTIGATION OF THE FLUTTER OF SEVERAL 
WINGS OF VARYING ASPECT RATIO, DENSITY, AND THICKNESS 
RATIO AT MACE NUMEERS FROM 0.60 TO 1.10 
By Raymond Herrera and Rotert E. Barnes 

SUMMARY 


An experimental Investigation has "been conducted to determine the 
flutter characteristics of several rectangular wings of variable aspect 
ratio at Mach ntmibers from 0.60 to 1.10. The wings were solid aluminum 
or steel and were 2, h, and 6 percent thlck- 

Flutter was encountered at both low and hi^ angles of attack. At 
low angles of attack, the flutter frequencies were of the order of the _ 
first natural bending frequencies e ven thcgugh the modes w ere of the / 
bending -torsion type. At hi^ angles of attack^ the mode^ were essen- 
tlally pure torsion and the frequencies were nearly equal to the first 
natural torsion frequencies. In most cases where stall flutter occurred, 
it occurred below Tnayimum static lift. 


INTRODUCTION 


Flutter may be classified according to the nature of the unsteady 
flow by which it is sustained. Accordingly, that class of flutter accom- 
panied by attached or potential flow has been termed classical flutter. 

It may occur in subsonic, transonic, or supersonic flow and may involve 
either single or multiple degrees of freedom. However, according to 
potential theory, classical flutter is not dependent upon angle of attack. 
A second type of flutter is associated with separated flow and is related 
to angle of attack. If flutter occurs at hi^ angles of attack it is 
often termed stall flutter. However, a separation-type flutter may occur 
at low or moderate angles of attack at transonic speeds and it is not 
always possible to establish a definite boundary between this type of 
flutter and classical flutter. 
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The status of research on classical flutter Indicates that there is 
a need for experimental research at transonic speeds . Separation-type 
flutter at high subsonic and transonic Mach manbers has became of interest 
because of its occurrence, on missiles (ref. l) and because of the need 
for data Trfiich are applicable to the design of high-speed aircraft having 
thin wings (ref. 2). The present Investigation was undertaken to provide 
such experimental flutter data for a series of rectsmgular wings of vary- 
ing aspect ratio. The tests were performed on the transonic bump of the 
Ames l6-foot high-speed wind tunnel at Idach numbers ranging from 0.60 to 
1.10, corresponding to a Reynolds number range from 1.7 to 2.0 million. 


SYMBOLS 


A 

El 

GJ 

la 


M 

V 

V 

bwct 

b 

g 

Sa 

Sh 

m 

^a 


^.g. 


'•e.a. 


aspect ratio, full span 
bending stiffness, lb-in. ^ 
torsional stiffness, Ib-in.^ 

polar moment of Inertia about the elastic axis per unit length, 
slug-ft^/ft 
Mach nimiber 
velocity, ft/sec 
reduced velocity 
wing semi chord 
damping coefficient 

structiiral dan^Jlng coefficient for the torsion mode 
structural damping coefficient for the bending mode 
mass per unit length, slugs /ft 

dimensionless radius of gyration of wing section, 

center of gravity location from leading edge, percent chord 

elastic axis location from leading edge, percent chord 


a 


angle of attack, deg 
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minimum angle of attack at 'wiiicli stall flutter occurred, deg 

mass density ratio, 

tunnel air density, slugs/cu ft 

air density at sea level, slugs/cu ft 

first natural torsional frequency, radlans/sec 

first natural "bending frequency, radians/sec 


MODEIS 


The seven rectangular vings tested had symmetrical NACA 64A-serles 
sections with thiclmess -to -chord ratios of 0.02, O.Olj-, or 0.06. The 
models had 0.5-foot chords and were constructed of solid steel or alumi- 
num. Two of the aluminum wings, percent and 6 percent thick, were 
modified hy cutting spanwise slots 0 . 005 - inch wide in the upper and lower 
surfaces to reduce the torsional stiffness. The models are designated 
in this report hy a number which Indicates the thickness (2, 4, or 6 per- 
cent) followed hy a letter to signify the material (a or S) . A letter X 
is used to designate a slotted wing. Thus, wing 4 ax indicates a 64 a 004 
aluminum wing with spanwise slots. 

The following physicstl characteristics are given for each wing in 
table I: mass density ratio, radius of gyration, bending stiffness, tor- 

sional stiffness, elastic axis, and center of gravity. "Values of "bojji 
and "bcija, based on the experimentally determined natural first bending and 
first torsional frequencies are given in table H for various aspect 
ratios of each wing. The structural damping coefficients of wings 2S, 

4S, 4a, and 4AX which were determined experimentp,lly are presented in 
table III. 


EXPERIMENTAL DETERMINATION OF THE ELASTIC AXES AND 
NATURAL FREQUENCIES 


The elastic axis of each wing was determined by noting the movement 
of a reflected ll^t beam while a wei^t near the tip was moved in the 
chordwlse direction. The distance of the weight from the leading edge 
for which torsional deflection of the wing could not be detected is given 
in table I for each model in percent chord. The wings were so mounted 
as to have their TnaxlTnum aspect ratios for these measurements. 
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The resonant frequencies were determined In the wind tunnel prior 
to testing. The equipment used for their determination is shown in the 
schematic and block diagrams of figure 1. The output of an audio oscil- 
lator was amplified and used to drive two 12 -inch speakers. The speakers, 
mounted above and below the model, were used to excite the model through 
two strings attached to the speaker cones. The output of a crystal vibra- 
tion pickup mounted at the wing tip and a reference signal from the 
driving oscillator were used to display Lissajous patterns on an oscil- 
loscope. These patterns permitted the operator to ascertain the resonant 
frequency. 

For the purpose of exciting the bending mode, the speakers were 
located in line with the elastic axis; for exciting torsion, the speakers 
were offset to each side of the elastic axis and their polarities were so 
related that a couple was produced. . 


EXPERIMENTAL DETERMINATION OF STRUCTORAL DAMPING 


The structural damping of wings 2S, 4A, and 4AX at various 

aspect ratios was determined by means of a free -oscillation technique. 
Bending oscillations were excited by the sudden release of a wei^t sus- 
pended from the tip at the elastic axis. A weight of 6.25 pounds was 
used on all wings except 2S for which a weight of 5 pounds was used. 
Torsional vibrations were excited by a sharp blow struck on the leading 
edge near the tip. This procedure unavoidably Introduced a small amount 
of bending which probably influenced the calculated damping coefficients . 

The motion of the wing tip was detected by a small crystal acceler- 
caneter whose mass was considered to have negligible effect on the results. 
The output of this accelerometer was recorded on an oscillograph. From 
these records the logarithmic decrement was determined for Intervals of 
about 20 cycles and the structural damping coefficient was calculated 
according to the following relationship (see ref. 3 )! 

g = i (logarithmic decrement) 

The coefficients so determined from several records were averaged and 
the results are presented in table HI. 


TUNNEL TEST APPARATUS AND INSTRUMENTATION 


The tests were conducted on the transonic bump of the Ames l6-foot 
high-speed wind tunnel. The wings were movable spanwise through clamps 
contoured to the wing profiles. Tie models were so positioned as to 
provide the desired aspect ratio and were then clamped by a force of 
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800 pounds olDtained tlirou^ the action of an air cylinder. The aspect 
ratio was varied from 2 to 6 for the and 6-percent -thick wings and from 
1 to k for the 2 -percent -thick wings. The maximum angle of attack was 
limited to 19 * 5 ° ky the angle -of -attack mechanism. 

The flutter frequencies were determined with the aid of a modified 
phonograph crystal cartridge mounted between the upper clamping block and 
the tunnel wall. The wing vibrations were transmitted to the pickup 
throu^ a short stiff wire resting on the upper surface of the wing. When 
the vibratione were at a constant frequency, the signal from the pickup 
was used in conjunction with a reference signal from an audio oscillator 
to display LissaJoTis patterns on an oscilloscope. 

A high-speed motion-picture camera (ilOO frames per second) was 
mounted outside the test section, opposite and slightly above the model 
and was used to photograph the upper surface and the wing tip when flut- 
ter occurred. The motion -picture records were used to aid in the quali- 
tative determination of the flutter modes. 


METHOD OF TESTING 


The Mach number range investigated was from O. 6 O to 1.10, correspond- 
ing to a Reynolds nximber range from 1-T to 2.0 million. The variation of 
velocity and relative density with Mach number is given in figure 2. At 
each test Mach number selected, the model was set at 0° angle of attack 
and so clamped as to provide its fflinlTmim aspect ratio. The angle of 
attack was then increased until (a) flutter was noted, (b) excessive sta- 
tic stresses or vibrations were encoiintered, or (c) a maxi mum angle of 
19*5° "was reached. The model was under continuous observation in order 
to prevent its failure and the observer made notes concerning its behav- 
ior. These notes are the basis for classifying flutter conditions as 
being mild or severe. When steady flutter occurred the frequency was 
detennined and in some instances a motion-pi ctixre record was taken. The 
angle of attack was then decreased to zero, the aspect ratio was increased 
hy 0.33^ and the procedure was repeated. Motion -picture records were 
taken of 21 cases of intennlttent flutter at low angles of attack; at 
high angles of attack, 25 steady stall-flutter conditions were recorded. 

The maximum allowable angles of attack were estimated on the basis 
of static wing loads and a yield strength of 70,CXX) Ib/sq in. for steel 
and 40,000 Ib/sq in. for aluminum. These limits were not rigidly adhered 
to and were subject to change during the test at the discretion of the 
observer. 

The effectiveness of the wing clamp was checked by retesting 
observed stall-flutter conditions (i.e., M, A, oj) while the wing clamp 
was rigidly bolted. The singles of attack at which stall flutter occurred 
with the wing clamp bolted were the same as those when the clamp was 
actuated by the air cylinder. 



6 


NACA EM A54A29 


RESULTS AMD DISCUSSION 

The results of the teats are presented in figures 3 to 9 for all 
wings tested. No flutter was obtained for aspect ratios other than those 
indicated in the figures. All flutter data presented represent either 
steady or intennittent flutter as indicated by a steady or Intermittent 
Llssajous pattern on the oscilloscope. 

In some Instances (fig. 4(a)) where mild flutter was noted at 
a = 0 °, the vibrations ceased with but a moderate increase in angle of 
attack and reappeared at higher angles as stall flutter. Where the flut- 
ter was more violent, particularly for the higher aspect ratios of wings 
2S and 2A (figs. 3 and 4), the vibrations persisted to higher angles of 
attack and limited the investigation to the angles shown. It was not 
possible to establish a definite boundary between classical and separation 
type flutter in all Instances. Only the steady flutter (figs. 3 to 9) at 
high angles of attack and the flutter which occurred initially at low 
angles ( 0 ° to 2 °) were considered to be subject to '.classification as 
stall flutter and low -angle -of -attack flutter, respectively. 


Flutter at Low Angles of Attack 


The determination of specific flutter points at low angles was not 
as decisive as -was the case at high angles. The motion pictures taken 
during the tests indicated that the flutter of wings 2A and 2S at small 
angles of attack -was of the coupled bending-torsion type in some Instances 
however, the predominant frequency was of the order of the first natural 
bending frequency. The vibrations of wing 2A (fig. 4) were quite severe 
from Mach number 0.80 to 1.10 for aspect ratios of 2.33 and 2.67. The 
vibrations were very violent and of large amplitude for aspect ratio 3*00 
and the predominant frequency from Mach nxmaber O .85 to 1.00 was approxi- 
mately 50 cycles per second. The model failed at M = 1,02, A = 3-00 
as the tunnel speed -was being reduced. The failure is believed to have 
been a result of structural fatigue rather than of divergent flutter since 
the model appeared to be more stable at Mach numbers above 1.00. 

The intermittent flutter noted for wing 2S (fig. 3) was near the 
first natural bending frequency for aspect ratios 2.67 to 3-33- The 
vibrations for aspect ratios 3-67 and 4.00 became progressively more 
severe with increasing Mach nimiber, so that the Mach number was limited 
to 0 . 94 . No flutter frequencies were determined for these aspect ratios. 

The vibrations of wing 4AX (fig. 7) at low angles were also near the 
first natural bending frequency for aspect ratios of 3*33^ 3-67^ and 4.00 
and did not appear to be in a coupled mode at any time. The vibrations 
of this wing when positioned at aspect ratios greater than 4.00 were too 
erratic to determine the frequency. 
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No lov -angle -of -attack flutter was encountered for wings 4S, 4A, 

6A, or 6AX. The intermittent flutter at moderate angles of attack which 
was noted for wing 4 a (fig- 6; A = 3-67# 4.00; a = 3-5°) and wing 6A 
(fig. 8; A = 6.00; a = 4,5°) was not considered to he of the low-angle- 
of -attack type since no flutter occurred at ot = 0°. 


Flutter at High Angles of Attack 


The stall flutter was of a steady torsional mode at or nesir the 
first natural torsional frequency of the wing, a phenomenon which sub- 
stantiates results of previous investigations. The predominant torsional 
mode was clearly evident in the motion pictures for various stall -flutter 
conditions. At a given Mach number, the intensity of the vibrations 
generally Increased with angle of attack to the limits of the test. How- 
ever, for wing 2S at A = 1.67 (fig- 3)^ it was possible to increase the 
angle of attack until a stable condition was again established. 

Several of the wings tested appeared to be less susceptible to stall 
flutter at the hi^er test Mach n\imbers. This trend was particularly 
evident in the case of wing 4AX (fig- 7j A = 2.00) for which steady flut- 
ter occurred at M = 0.90 (a = 16 . 5 °) ^ intermittent flutter occurred at 
M = 0.94 (a = 15 . 5 ° bo 19-5°)j no flutter was noted from M = O .98 

to M = 1.10 at angles of attack up to 19-5°- A similar but less pro- 
nounced trend was noted for this wing at A = 2.33 and 2.67, as well as 
for wing 6AX (fig. 9; A = 3-33 to 4.67) # and wing 2S (fig. 3; A = I. 67 ). 

As mentioned previously, the occurrence of steady stall flutter was 
pronounced and easily determined by the Lissajous patterns on the oscil- 
loscope. The minimum angle of attack at 'vdiich steady flutter occurred 
(of ) has, therefore, been plotted as a function of Mach number in figures 
10(a) to 15 (a) to summarize the effects of Mach nuniber and aspect ratio 
for all wings except 6 A. This wing exhibited only a mild form of inter- 
mittent flutter within the angle -of -attack range investigated. 

Data from part (a) of figures 10 to 15 are presented in a different 
form in part (b) of these figures to Introduce the conventional flutter 
parameter v/btoo, s-s a function of ct^- The parameter V/btuico decreased 
rapidly with increasing angle of attack. A value of V/btOa = 1 has been 
suggested in reference 2 as being useful in estimating the minimum stall- 
flutter velocity of thin wings. Explicit verification of this criterion 
is not possible, however, since no data were taken at Mach numbers less 
than 0.60 or at angles greater than 19 - 5 ° • 

In figure 16 , oj* is plotted as a function of aspect ratio at 
M = 0.60 to indicate the effect of a reduction in torsional stiffness 
on Of for the 4-percent and 6 -percent aluminum wings. As noted previ- 
ously, wing 6 A exhibited only mild intermittent flutter up to the Tnaviimmi 
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angles of attack aHowatle from mechanical or strength considerations. 
These values of a, shovn in figure l6 for this wing, are considered to 
be conservative estimates of the true values of The data as pre- 

sented indicate that reduction of torsional stiffness had a more pro- 
nounced effect for the 6-percent wing than for the 4- percent wing. A 
similar comparison at other Mach numbers is not possible, however, because 
of the lack of data for the solid wings. Wing 4A failed during the test 
before sufficient stall-flutter data were obtained. The data for wing 
6 a at Mach numbers higher than 0.60 were not deemed sufficiently reliable 
to continue the comparison. 

Stall flutter has been associated in the past with the static lift 
curve (see ref. 4), This association led the authors of reference 4 to 
suggest three possible causes of stall flutter; (l) static instability 
due to the negative slope of the lift curve at angles beyond the stall, 

(2) dynamic instability due to hysteresis at the stall, and (3) insta- 
bility resulting from Karman vortex excitation. 

In figures 17, l8, and 19 the suigles of attack at which flutter was 
observed have been noted on static lift curves taken from reference 5» 

These static data were obtained from tests conducted in the same facility 
and on models of the same size as those of the subject investigation. It 
must be noted, however, that these ctirves are for wings having symmetrical 
MCA 63A-series sections rather than the 64A-series used in the present 
investigation. Comparison of available 63-series and 64-series -airfoil- 
section data in references 6, 7, 8 indicates that the static data 

used are Indicative of the true static characteristics. 

The limited number of available data as shown in figures 17 to 19 
indicate that with respect to the three possible causes of stall flutter 
cited above, only the steady flutter of wing 4AX at the lower aspect ratios 
(fig. l8) and the intermittent flutter of wing 6A, A = 6.00 (fig. 19) may 
be attributed to static instability due to the negative slope of the lift 
curve. The remaining cases of steauiy flutter can then be considered as 
probably being a result of the hysteresis loop described by the lift curve 
since the stall -flutter frequencies were considerably lower than the vor- 
tex frequencies calculated by the method given in reference 4. 


CONCLUDING REMARKS 


Flutter at low angles of attack occurred only for the 2 -percent- 
thick steel and alumintmi wings and for a 4 -percent -thick aluminum wing 
slotted spanwise to reduce the torsional stiffness. TSie predominant flut- 
ter frequency was of the order of the first natural bending frequency. 

Stall flutter occvarred at frequencies -which were in good agreement 
with the experimentally determined first natural torsional frequencies. 


iliWPIIlFPII HU I 
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The angles of attack at -which stall flutter occurred have heen 
indicated on the static lift curves of a series of rectangular -Hrings 
having symmetrical NACA 63A-series sections.- In the majority of cases, 
flutter occurred at angles helow those for TnavtmiTm lift. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aesronautics 
Moffett Field, Calif., Jan. 29 f 195^ 
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Wing 

section 

Material 

Mass 

density 

ratio, 

ti 

Dimensionless 
radius of 
gyration, 

CL 

Bending 

stiffness, 

El, 

Ib-in.^ 

Torsional 
stiffness, 
GJ, 
lb -in 

^.a.^ 

per- 

cent 

chord 

per- 

cent 

chord 

23 

64A002 

Steel 

108 

0.460 

1 . 15 x 10 ^ 

1.84x10* 

37.3 

43.5 

2A 

64AOOe 

Aluminum 

38 

A99 

0.403 

0.594 

38.3 

43.5 

Its 

61tA004 

Steel 

216 

.410 

9.24 

14.77 

4o.o 

■ 43.5 

4a 

61tA004 

Aluminum 

76 

.409 

3.23 

4.77 

43.3 

43.5 

^4AX 

64A004 

Aluminum 

76 

.571 

3.23 

3.51 

43.5 

43.5 

6a 

64a006 

Aluminim 

ll4 

.425 

10.92 

16.1 

46.0 

43.5 

®6AX 

64aoo6 

Aluminum 

ll4 

.547 

10,92 

10.6 

40.5 

43.5 


^Spanwlse Blots having a depth of 55 percent of the ■wlng-eectlon ordinate located 

at chordwlae Intervals of 0.3 inch between 5- 90-percent chord. 

^Spemwiae slots having a depth of 35 percent of the wing-section ordinate located 

at intervals of 0.3 inch between 5- and 90 -percent chord. 

Note: Values of mass density ratio n are based on the density of air at sea 
level. For variation of test density with Mach number see figure 2. 
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TABLE II.- STILL-AIR FOECED-VIBRATIOH CEARACTERISTICS OF WIHQS TESTED 


I 



Aspect ratio 

Wing 

2 

3 

4 

5 

6 


isa 



IS 

^»2hOn 





“b“a 

2B 

U7-8 

323.5 

55.2 

213.5 

34.1 

158.5 

. - » 


^ « M 

. _ . 

2A 

119.5 

311.0 

56.6 

207.0 

35.3 

152.4 

~ - 

*■ PA H 

•A PA *■ 

AA i« A* 

hs 

M» M M 

PA n M 

94.2 

436.0 

59.8 

342.0 

4l.0 

278.0 

29.2 

231.0 

kk 

•V w w 

. - - 

92.6 

432.0 

55.8 


38.0 

262,0 

28.8 


ItAX 

263.0 

575.0 

116.8 

383.0 

66.7 

284,0 

48.2 


37.2 

195.0 

6A 

- - - 


152. 

644.0 

88.8 

474.0 

60.0 

400.0 

42.6 

330.0 

6 AX 

A*. «■ M 

■ 

i57.a 

534.0 

91.1 

393.0 

60.0 

311.0 

45.4 

259.0 


ving Bemlchord, 0.25 ft. 
first natural 'bending frequency; radians/sec. 
®0i\x first natural torsional frequency, radlans/Bec. 




mBIE III.- BPERmEAUT EBIBKQIED BraKITOAI, BAMPim CQEFVICIEinS 



Aspect ratio 

Wing 

2 

3 

4 

5 

6 

\ 

% 

\ 

X 

"8h 

X 

X 

X 

X 

X 

2S 

4s 

4 a 

4 ax 

1 

1 

0.007 

.013 

.007 

0.004 

.004 

0.005 

.005 

.007 

.006 

0.004 

.003 

.004 

0.006 

.006 

1 

0.005 

.007 

till 

1111 

lilt 


structural danqiing coefficient for the hendlng mode, 


structural damping coefficient for the torsion mode. 
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Oscilloscope 


Audio 

oscillator 


Audio 

amplifier 



Figure 1.- Schematic and block diagrams of apparatus used for the 
determination of the natural frequencies. 
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(b) A = 3 . 00 , 3.33, 3.67, ^.00 
Figure 3.- Concluded. 








(a) A » 1.67, 2.00, 2.33, 2.6? 
Figure 4.- Flutter obaervatlonB, wing 2A. 
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(b) A = 3.00 
Figure 4.- Concluded. 
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(b) A »= 5.00, 5 
Figure 5*- Co: 
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O lnf9railtfot flutter 

? Sf9atfy flutt9r 
Maximum a t99f9<l 
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(b) A - 4.33, 4.67 

Figure 6.- Concluded. 
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(a) A « 2.00, 2.33, 2.67, 3-00 
Figure 7'" Flutter observationa, vlng 4AX 
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Figure 10.- Effect of aspect ratio and Mach number on the stall flutter 

of wing 2S. 
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Figure I?.- Stall flutter of the 64A002 wings related to the static lift 

curves of a 63A002 wing. 
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Figure iB.- Stall flutter of the 64AOOh wings related to the static lift curves of a 63A004 wing. 
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Figure 19.- Stall flutter of the 64AD06 vlngs related to the static lift curves of a 63AD06 wing. 








